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Abstract 

Exact relations between gauge-invariant vacuum correlators in QCD 
are derived. Derivatives of the correlators are expressed in terms of 
higher orders correlators. The behaviour of the correlators at large 
and small distances due to these relations is discussed. 
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Understanding and description of nonperturbative interactions is one of the 
most important problem in the modern held theory. Some progress on this 
way has been achieved with the method of vacuum correlators (MVC) [1]. 
Hi. Nonperturbative interactions are expressed in this method through 
gauge-invariant correlators of the gluon held strength operators, acting at the 
different space-time points. All gauge-invariant quantities like, for example, 
string tension and hadron masses may be written in terms of these correlators 
[§] . On the other hand, the latter can be computed on the lattice |4j , and thus 
play a role in the theory as intermediate step for representing the dynamical 
information. 

To have a self-consistent theoretical scheme for MVC, it is necessary to 
derive the equations for correlators, or, at least, some relations between dif¬ 
ferent correlators. 

Concerning the hrst, a set of exact equations has been proposed [|J, which 
obtained by the stochastic quantization method, but, unfortunately they look 
quite difficult to be solved and there are still no physical results derived in 
this framework. 

We choose the second way in this paper and derive simple relations, yield¬ 
ing some information about the behaviour of the correlators at small and large 
distances. In the previous paper the equation for 2-point and 3-point cor¬ 
relators was found. The analysis of this equation will be continued here and 
the higher order correlators - 4-point ones are also included. 

We are going to study gauge-invariant Green functions ( n-point corre¬ 
lators) of the following form: 

=< Tr{F^ lVl (z 1 )^z^z 2 )F^ V2 {z 2 )...F^ nVn {z n )^(z n ,z 1 )) > ( 1 ) 

here F fiU (z) = F^ u (z)t a - is gluon held strength and nonabelian phase factors 
(parallel transporters) <h(^, z') = P exp i / 2 2 A^dx^ - are crucial for gauge- 
invariance of the correlators. Vacuum averages are defined as follows: 

< 0(A) >= J 0(A) 

Any n-point correlator may be constructed using invariant tensors £^a P , 
vectors (zi — Zj) /1 and a set of scalar functions, depending on relative coor¬ 
dinates. For the simplest nontrivial 2-point correlator one has two scalar 
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functions |T|: 


Agli,u 2 , 2 =< Tr{F^ 1 (z 1 )^(z 1 ,z 2 )F tl2U2 (z2)^{z 2 ,z 1 )) >= 


2 (^/U2^1^2 Z U2 5 uip2 ) T i^U2^H1H2 ^^2 ^ 1 ^ 2 )^ -^l(~l ^ 2 )“^ 


T 2 ^l 1 l l/ l a b^'^2^2^ D(yAj 1 ^ 2 ) 

Function -D(^) as well as ^(z) may be written as a sum of two parts: 


( 2 ) 


D(z) = D pert (z) + D np (z ) (3) 

The first part D pert (z) —> cx) if —> 0 and physically represents the exchange 
of perturbative gluons. Notice here, that one-gluon exchange does not get 
contribution to the D pert (z), the contribution to Di pert (z) is of the following 
form: 

^ . . 16a s (z 2 ) . 9 . 

Dipert{z) = + 0(cc;). (4) 

At the same time for the theories with a nontrivial vacuum like QCD, there 
is the second part D np (z), which has no perturbative interpretation. Lat¬ 
tice simulations show [||], that the function D np (z ) for |;z| > 0.2 fm falls 
exponentially: 

D np (z) ~ e~^a 

with the correlation length T g = 0.2 fm, the behaviour of D np (z ) for smaller 
z is still unknown. At the point z — 0 nonperturbative components are 
expressed through the gluon condensate ||: 

^\i}ij\,ij.2V2 (^) =< ^ 2 ^ 2 ) > = 

= 2 liiviabtfi 2 ^ 2 db ' ( D(0) + £>i(0) ) (5) 
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Defined as in (2), z’) depends on a contour, connecting the points z and 

z'. It is convenient to consider a set of correlators, where transporters $ 
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connect neighbouring points along the choosen trajectories. We will consider 
correlators with the points, lying on a fixed straight line, this choice is natural 
for 2-point correlator and we follow it for higher correlators too. 

The fields under consideration are nonabelian, that’s why the derivative 
of the correlator with respect to z a is not equal to the correlator of D a F^ v (z) 
with another fields. The difference arises from the differentiation of the 
contour and it is crucial for the nontriviality of the picture. Nonabelian 
transporters 

f z ' 

3>(z, z') = P exp i / A^dx^ (6) 


are differentiated according to the rule 
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d$(z, z') 

dz^ 


i$(z, z^A^z') + i(z' - z) p I fn (z, z') 


( 7 ) 


where notation is used: 

I py (z,z')= / da a $(^, z + a{z' — z))F fry (z + a[z! — z))- 
Jo 

■ <f>(z + a(z' — z), z') (8) 

Analogously 


~~ = z') + i(z' - z) p I fn (z, z') (9) 

I fry (z,z')= / da ■ a^(z, z 1 + a(z — z'))F p r Y (z'+ a(z — z 1 ))- 
J 0 

• $>(z' + a(z — z'), z 1 ) (10) 

We shall often omit the arguments of transporters below for the sake of 
simplicity . Let us consider the derivative of the 2-point correlator (2), using 
(7) and (9) one has: 

=< TriF^zjQDtF^fo)*) > + 

OZ 2 £ 

+i(z 2 - Zi) a (^< Tr(F lilvl (z 1 )I a z(z l ,Z 2 )F fi2U2 (z 2 )®(z 2 ,z l )) > - 
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- < Tr{F^ Vl {z 1 )^{z 1 ,z 2 )F IJl2V2 {z2)I a ^{z2,z 1 )) > 

The l.h.s. of this equations is: 


£ P2V 2 


a (2) 


cb 


2? 


cb 2 




( 11 ) 


( 12 ) 


In such a manner one obtaines a connection between the 2-point and 3-point 
correlators: 


-mvKTp 


CW ^ = ^ £ m2^p(< rr(F M1I/1 (^ 1 )/ (T? (^ 1 ,^ 2 )F M2 ^(z 2 )<I)(2:2,-i)) > 


- < Tr(F lxl „ 1 (z 1 )^(z 1 , z 2 )F IJi 2 U 2 (z 2 )Ia^(z 2 , > j (13) 

We have taken into accout the Bianchi identity £ P 2 U 2 ^ p D^F^ 2 V 2 {z) = 0 and 
denoted z 2 — z\ — z. 

This relation is true for the perturbative as well as for nonperturbative 
parts of the correlators. We will focus our attention on the latter, assuming 
regular behaviour if their arguments tend to zero. The expression (13) is 
simplified if z —> 0, because in this case I a ^(zi,zi) = I a ^(zi,zi) = ^F a ^(zi). 
Thus one obtains: 


dD{z) 
dz 2 


2=0 


/ 


abc 


96 


< F a F b F c > 

pi VI V1V2 V2P1 


(14) 


Relations (13) and (14) have been derived in [|6|. Let us take up a question, 
what tensor structures of the 3-point correlator contribute to the r.h.s. of 
(14). There are two independent Croneker tensors in this case, namely: 


A 0) 

P1V1,P2V2,P3V3 


< Tr ( F p , <D> F /12 „ 2 4> F /J3 V3 $) >= 


g £piviab^p 2 v 2 bc£paused ' D3 T 

l ( \ 

g I 1 ~piviv 2 a,£/i 3 v 3 p 2 a £pivi^a^P 3 V 3 V 2 C 1 T permut. J 

+ nonkroneker terms. 


• D2 + 


(15) 


It is easy to see, that the part, proportional to D s does not contribute, while 
D 2 does: 


rabc < pa 7 7 b pc 

J ' piv 1 ± viv 2 v 2 pi 


>= 48T> 2 (0). 


(16) 
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or 


dD(z) 


dz 2 


D 2 (0) 


(17) 


z =0 


Formulas (16), (17) have a simple interpretation in terms of the dual 
Meissner effect. Introducing three-dimensional notations for chromoelectric 
and chromomagnetic fields F= £ a( g 7 iL“ , left hand side of (16) 
may be presented as: 


r 


abc 


< F^F^F^ >= r%3 7 (3 < E*E%H$ > + < >) (18) 


Thus a nonzero vacuum averages of triple correlators imply, that spontaneous 
magnetic fluxes may be created due to electric or magnetic lines splitting. 
The existence of noncorrelated magnetic fluxes of this type is equivalent to 
the presence of monopole currents and is the key point for the possibility of 
confinement in such vacuum. Notice also, that colour structure of the 3-point 
condensate is determined by antisymmetric constants f abc , symmetric colour 
term is absent: 


d 


abc 


< F a F b F c >= 0 

Ml*'! M2 ^2 M3 ^3 


or 


fabc 

< pa pb pc >= £ _ A (3) 

M 1^1 M2 V2 M3 F3 g Ml^l ? M2^2,M3^3 


( 0 ) 


The equation (13) will be examined here in more details. One can prove, 
that it is equivalent to : 


dD{z) 

dz 2 


1 

2 


tdt (^D 2 (z — t,t) + D 2 (-t,t — z) S j 


(19) 


Let us make an assumption, that function D and D 2 may be parametrized, 
taking into account their cyclic symmetry, as: 

D(z) = D(0) h{z)h{—z)\ D 2 (z — t,t) = D 2 (0) h 2 (z - t)h 2 (t)h 2 (—z) 


Then, it is easy to check, that the solutions of the following equation satisfy 
also (19): 


h{-z) 


dh{z) 

dz 
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dt h 2 (t ) • h 2 {z 


t) 


( 20 ) 
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where a 2 = — 2 olo)' > • Naively one can think, that h(z ) = h 2 (z) for all z. 
Equation (20) may be solved in this case and the regular at z = 0 solution 
has the form: 

h(z) = 2 

az 

where Ji(oz) - Bessel function of the first type. This behaviour contradicts 
the lattice measurnrents [|4|, it means, that the relation h(z) = h 2 (z) does 
not hold in the real gluodynamic vacuum. Still this solution leads to the 
def ini te value of the string tension, a: 


a = 


J d 2 z D(z) 

o 


M-P(Q )) 2 

d 2 ( o) 


The important point here is the negative sign of the fraction he. a 

should be real, overwise correlation functions would unphysically increase 
with a distance. Indeed, there are several grounds to think, (see 0 and ref. 
in |Q]) that the 3-point condensate is negative in the real world. 

Let us consider a more realistic case h(z) ^ h 2 (z) . Lattice results allow 
to look for the solution in the form: 

1*1 °° d ( z \ 11 
D(z) = D( 0)e-^ f 

n=o n - \ ± aJ 


M 


D 2 (x,z - x) = D 2 ( 0) e T a £ 


X 


l*-*l 


771=0 


m\ l T n 


■ e 


£ 

m =0 


b m (z- x 


ml 


T 

-L n 


( 21 ) 


The substitution of (21) in Eq. (19) gives a relation between d n and b m \ 

(gT g y 

2 m= 0 


dn-\- 2 — 1 


n+1 


^ ^ bm bn—r 


( 22 ) 


with the obvious conditions do — bo — d\ — 1. There are infinitely many 
sets of coefficients b m ,d m , which solve the above equation. What solution is 
realized in the real world is a question for future investigations. 
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It remains to be seen whether the above picture is compatible with the Gaus¬ 
sian model for QCD vacuum |Q], ||. In the latter all correlators of the odd 
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orders, including A® are assumed to be zero, while the even correlators 
are factorized in a product of the 2-point ones. Among main arguments in 
favour of it is the fact, that potential of confinement for higher representa¬ 
tions of a gauge group is proportional to the quadratic Casimir operator in 
the Gaussian vacuum. This fact is confirmed with a good accuracy by lattice 
computations (see discussion in 0). The Gaussian model for QCD vacuum 
was introduced in [11] and intensively used in [[?J for the high energy scattering 
amplitudes. Eq. (19) shows clearly, that such picture is not selfconsistent 
for any nontrivial function D(z). Thus one should introduce an extended 
Gaussian ensemble and postulate factorization in a product of the 2-point 
and 3-point correlators. If 3-point correlator is zero, this definition coincides 
with the ordinary Gaussian ensemble, because all odd order correlators are 
also zero, in the opposite case a new stochastic ensemble arises, it is natural 
to call it the minimally extended Gaussian ensemble. At the same time even 
in the ordinary Gaussian vacuum function D(z) may be a nontrivial one. We 
shall deal in this case with the second derivative of the 2-point correlator at 


z = 0: 


0 2 A ( 2 ) 


£ Hivipn^ /i2^2?7 ^ 


X — 

%(J %<f) 


2 < Tr{F llll/1 [Fp^F IJi2l , 2 \) > + ^ < Tr{F lilVl [D p F (T ^F^ 2l , 2 ]) > 

^6 < Tr(F lllul F^p[F a £ F^ 2 „ 2 ]) > +6 < Tr(F M1I/1 [F M2!/2 F CTe ]F 0p ) > + 


+ < TriF^F^F^F^) > + < >) 

Taking into account the following: 


(23) 


d 2 AW(z) 

dz 2i dz lp 


d 2 D(z) 

(< dz 2 ) 


( Zq 


5-y V z 2 ) - 24 ——— 5, 


dD 
d? 


(24) 


one has: 


-48F"| 0 = lim 


tz 

^ 2 ^ 1 ^ pv e < Tr(F^ lUl [D p F a £, Fp^]) > 


(b <Tr(F^ lvl [F^ ]) > +6 < Tr{Fp lV1 [Fp 2V2 F (7i :}F 4>p ) > + 












+ < Tr(F Mm [F^i^ p ]F M2l/2 ) > + < Tr(F mui 

Fp 2 f2 [F(j>pF a (:]) > 


(25) 


The first and the second terms in the r.h.s. of (23) are absent in the 
Gaussian vacuum, 4-point condensate has to be introduced for the analysis 
of the others: 


A< 4 > =< TriF^F^F^F^) >= D~ ct + 


+ («0 ' & p\v\ab ^ p,2V2bc ^^ 3^3 cd T Oi 1 ' £ p\v\ab E-psi'sba £p2V2cd & /^i^dcT 


+ «2 • £ M i M 2 afe ba e^cd £ p 3 v 3 dc) • D c [ auss + nonkron. terms (26) 

here the connected 4-point function D^ onnect should also be zero in the Gaus¬ 
sian vacuum. 

Using (5): 


s :ab 


^ Fp\V\Fp 2 u2 ^ O ^TlI / lPfUp 2 P 2 PC r (D(0) + DM) 


(27) 


it is easy to get: 


«. = <* = * a , = a*-. = (« 


(28) 


Straightforward calculations lead to the result: 


d 2 D(z ) 


(dz 


2\2 


z =0 


7 ( D^ + D^Oy 
2 1 4 


(29) 


We see, that even for the Gaussian ensemble the presence of the nonper- 
turbative condensate 

D(0) + Di(0) ^ 0 (30) 

automatically provides the dynamical dependence on z for function D(z). 
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The relation (29) demonstrates the existence of a scale, which determines the 
behaviour of A^ 2 ^(z) at small distances and it may be quite different from 
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the asymptotic correlation length T g . Indeed, taking from |4|] D\(z) ~ \D(z) 
and defining D( 0) + Di(0) from the standard gluon condensate [|] , we get 

hkS ss 0.012 Gen 1 

B(0) 

whence it follows that (parametrizing function D(z) near zero as a gaussian) 

D(z) = D( 0)e-^ 2 / r o 

and r 0 ~ 0.6 fm. This value of r 0 is larger than T g « 0.2/m, characterizing 
the exponential asymptotic decreasing of the correlator for z > 0.2/m: 

D(z) = D{ 0)e~ WT ° 


if z > 0.2 . 

This situation is possible in the so called string limit: 


T g (D{ 0) + -D^O)) ~ cr = const 


with T g —> 0. It is seen from (29), that also: 


Tt d 2 D(z) 
D(0) + D!( 0) (dz 2 ) 2 


0 


z =0 


(31) 


Substituting in (31) values for D( 0) + T ) 1 (0) and T g from [Q], we get on the 

l.h.s. 0.04 ~ ^ , that makes string limit in the nonperturbative QCD 

credible. 

Let us make another important notice about the Gaussian approximation 
for the gluodynamic vacuum. The Taylore expansion for even function D(z) 
is as follows: 


s rc(z ) 

h (* 2 )" 


z=0 


2 n 


The coefficients of the above expansion are expressed through the (2 + n )- 
th and less orders condensates. We have seen that for n — 1 and n — 2 
the derivatives at z = 0 are zero and positive respectively for the Gaussian 
ensemble. It may be interesting to check, are there negative coefficients in 
this Taylor expansion for the Gaussian vacuum or not. The latter would be 
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obviously incompatible with the observed decreasing of the function D(z) 
with z and would be a rigorous proof of nongaussian nature of gluodynamic 
vacuum. 

The behaviour of the function D(z ) for z < 0.2/m is unknown at the 
moment. Its determination and higher correlators measurments would take 
a chance to check the relations (13), (14), (22) and make a conclusion, what 
type of vacuum is realized in gluodynamics and whether one can approximate 
this vacuum as the Gaussian or the minimally extended Gaussian ensemble 
of the stochastic fields. We have said nothing about perturbative parts of cor¬ 
relators and their connection with the nonperturbative ones. This important 
question will be considered in its own right. 
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